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Abstract 

1 the whole complex 

is studied, assuming that 



Strong asymptotics on the whole complex plane of a sequence of monic Jacobi polynomials p^"'"''") 



lim ^ = .1 , lim ^ = B , 

with A and B satisfying A > —1, B > —1, A + B < —1. The asymptotic analysis is based on the 
non-Hermitian orthogonality of these polynomials, and uses the Deift/Zhou steepest descent analysis 
for matrix Riemann-Hilbert problems. As a corollary, asymptotic zero behavior is derived. We show 
that in a generic case the zeros distribute on the set of critical trajectories F of a certain quadratic 
differential according to the equilibrium measure on F in an external field. However, when either a„, j3„ 
or a„ + Pn are geometrically close to Z, part of the zeros accumulate along a different trajectory of the 
same quadratic differential. 



1 Introduction 

We consider Jacobi polynomials p^'^"'^") with varying negative parameters An and Bn such that 

-l<yl<0, -1<B<0, -2<A + B<-1. (1.1) 

We will obtain strong asymptotics as n — )■ oo of p^^^"'^") miiformly for z in any region of the complex 
plane and uniformly for A and B in compact subsets of the set of parameter values satisfying (11.1(1 . Since the 
asymptotics is uniform in A and B, we also find the asymptotics for general sequences of Jacobi polynomials 
p{a-n,Pn) g^(,j^ t\ie.t the limits 

A= lim — and S ^ lim ^ (1.2) 

n — >oo fi n — >oo n 

exist, and satisfy (|l.l|l . From the asymptotics of the polynomials we will also be able to describe the limiting 
behavior of the zeros. 

From the point of view of behavior of zeros, the Jacobi polynomials with varying parameters q;„, /3„ 
such that 1(1.2(1 and 1(1.1(1 hold arc the most interesting general case. Indeed, Martmez-Finkelshtein et al. 
PK] distinguish five cases depending on the values of the limits (11.2(1 (cf. Fig. The first case is the case 
where A, B > 0, which corresponds to classical Jacobi polynomials with varying positive parameters. These 
polynomials are orthogonal on the interval [—1,1], and as a result all their zeros are simple and belong to 
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the interval (—1,1). The asymptotic behavior of Jacobi polynomials with varying positive parameters is 
discussed in j^I El ^1 ^1 IZHl IZE]- We also consider the parameter combinations B > 0, A + B < —2 and 
A > 0, A + B < —2 as classical. Indeed, the transformation formula 

see [321 §4 .22], expresses a Jacobi polynomial with parameters a and /3 satisfying a + f3 < —2n and /3 > — 1 
directly in terms of a Jacobi polynomial with positive parameters. It follows that 1)1. 3|l reduces the study 
of Jacobi polynomials with varying parameters a„ and /3„ such that the limits 1)1. 2|l hold with B > and 
A -\- B < —2 to the study of Jacobi polynomials with varying positive parameters. The analogous formula 

shows similarly how to reduce the case A > Q and A + £? < — 2 to the classical case. 




Figure 1: The five different cases in the classification of Jacobi polynomials with varying parameters according 

to ESJ. 

The second case in the classification of corresponds to limits A and B in 11.2|l satisfying one of the 
three combinations A < —1, A + B > —1, or B < —1, A + B > —1, ot A < —1, B < ~1. In this case the 
zeros accumulate along an open arc in the complex plane. Their asymptotic distribution was found in |25| 
in terms of the equilibrium measure in an external field (cf . |29| ) . The approach followed there was based on 
the non-hcrmitian orthogonality of the Jacobi polynomials with these parameters. See |2()) for an overview 
of non-hcrmitian orthogonality properties of Jacobi polynomials with general parameters. 

The remaining cases correspond to combinations of A and B values such that one or more of the inequal- 
ities — 1 < A < 0, —1 < B < 0, and —2 < A + B < —1 are satisfied. In these cases, the zero behavior is 
more involved due to the possible occurrence of multiple zeros (at ±1 only) or a zero at oo (which means a 
degree reduction). To be precise, if a = —k is a negative integer with k £ {1, . . . , n}, then we have (see |32l 
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formula (4.22.2)]), 

so that Pn ^'^'^ has a zero at 1 of multipHcity k. Similarly, li fi = —I with I e {1, . . . , n} then Pi"'~'^ has 
a zero at —1 of multiplicity I. A degree reduction may occur when a + /3 is a negative integer, namely if 
a + P = —n — fc — 1 with k G {0, . . . , n — 1}, then 

^" " r(fc + a + 1) n! ^'^ ^^-^^ 

see 122 Eq. (4.22.3)]; see §4.22 of 32 for a more detailed discussion. Now assume we have varying parameters 
Oin, Pn such that the limits (|1.2() exist. If — 1 < A < 0, and if the q;„ are integers, then we have for each n 
large enough, that p^^'f^"^ has a multiple zero at 1. In the weak limit of the zero counting measures this 
corresponds to a point mass |^| at 1. Similarly, if —1 < _B < 0, and if the /?„ are integers, then we have in 
the limit a point mass \B\ at —1. Finally, if —2 < A + B < —1 and a„ + /3„ are integers, then we have in 
the limit a point mass 2 + A + B at infinity. 

The classification of the remaining cases in depends on the number of inequalities — 1 < A < 0, 
— 1 < i? < 0, —2 < A + B < —1 that are satisfied. The third, fourth and fifth case correspond to 
combination of parameters A and B such that exactly one, exactly two, or exactly three, respectively, of the 
inequalities are satisfied (cf. Fig. In these three cases the limiting behavior of zeros will be very sensitive 
to the proximity of a„ (if -1 < A < 0), /3„ (if -1 < i? < 0) or a„ + /?„ (if -2 < ^ + S < -1) to integer 
values. For Laguerre polynomials the same phenomenon was analyzed recently in |21| . 

Since all three kinds of singular behavior can occur in the fifth case, this is the most interesting case and 
that is the reason why we consider it here. The other cases can also be treated with our methods. Fig. ^ 
shows the behavior of zeros which is typical for case 5. From the figure it appears that the zeros accumulate 
on a contour consisting of three analytic arcs. From our analysis below it follows that this is indeed the 
case, provided that the parameters are not too close to integers. We identify the curves as trajectories of a 
quadratic differential as well as the limiting density of the zeros on the curves, see Theorems 12.31 and 12 .41 for 
the exact statement. To be able to explain the remarkable zero behavior was the main motivation for the 
present work. 

We remark that the different possibilities within the cases 1, 2, 3, and 4 can be transformed to one 
another using the transformation formulas H1.3I) . (|1.4|) for Jacobi polynomials. It is interesting to note that 
case 5 is invariant under these transformations, see ,25, . 

We also remark that the transitions between the five cases (i.e., A — 0, A — —1, or B = 0, B = —1, 
A + B = —1 or A + B = —2) will present additional difficulties. These are the non-general cases, in contrast 
to what we call the general cases 1-5. The zero distribution in some of these cases has been studied by Driver, 
Duren and collaborators (see also a recent survey [33] on the large parameter cases of the hypergeometric 
function). In the case p^'^^+i'"""!) ^ k £ N, has been analyzed, corresponding to A = k € N and 
B = —1; this result was generalized in |15| using a saddle-point method to allow k to be any positive real 
number. Case pi""*"^' " has been studied in ^^l- In general, these works establish the accumulation 
set of the zeros but not the limiting distribution. Trajectories of the zeros of the Gegenbauer polynomials 
p(-n-b,-n-b) ^j^j^ fixed n ss b varies from —1/2 to — oo have been described in HJ. 

The rest of the paper is organized as follows. The main results are stated in Section 2. We start defining 
the basic configuration on the plane used in the description of the zero (Subsection 2.2) and strong (Subsec- 
tions 2.3-2.4) asymptotics of the polynomials. In Section 3 we prove two technical lemmas. The cornerstone 
of our approach is the matrix Ricmann-Hilbert problem formulated in Section 4; the transformations of this 
problem in the framework of the Deift-Zhou steepest descent analysis (Section 5) are used in Section 6 to 
prove the main results of the paper. 
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Figure 2: Zeros of -Pioo ^^^^ — ""^^ +1*^ P — ~80 +10 '"^j which corresponds to Case 5 in Fig.^ 

The Deift-Zhou steepest descent method for asymptotics of Riemann-Hilbert problems was introduced 
in [H] and apphed first to orthogonal polynomials in |H], see also (Hj- We use an adaptation of the 
method to orthogonality on curves in the complex plane. The optimal curves are trajectories of a quadratic 
differential and they were used for steepest descent analysis of Riemann-Hilbert problems first in |2| and 
later in H [T3 EH 1^ El ■ 



2 Statement of results 

2.1 Geometry of the problem 

We assume A and B satisfy the inequalities Hl.l|l and define 



B^-A^±At^{A+l)iB + l){-A-B-l) 

iA + B + 2)^ ■ ^ ' 

Because of the inequalities Hl.lll we have that all factors in the square root in 12.1(1 are positive, so that 
C+ G C"*" = {z &C : Imz > 0} and C- is the complex conjugate of C,+ . 

Regardless of the branch of the square root and of the path of integration we choose, the set 

r ^ r(^'^) := (z e C : Re ~ "^^f ~ dt = o] (2.2) 



is well defined, and consists of the union of the critical trajectories of the quadratic differential (cf. PJj) 



Lemma 2.1 We have that T is the union of three analytic arcs, which we denote by T^, Tc, and Tj^. All 
three arcs connect the two points and intersect the real line in exactly one point, in such a way that each 
of the intervals (— oo,— 1), (—1,1), (l,oo) is cut by one of the arcs. 
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The contour T is oriented as indicated in Fig.|21 That is, and Tc are oriented from to C_, and Tji 
is oriented from C_ to The orientation of F induces a + and — side in a neighborhood of the contour, 
where the + side is on the left while traversing F according to its orientation and the — side is on the right. 
We say that a function / on C \ F has a boundary value for i G F \ C_} if the non-tangential limit 

of f{z) as 2 — * t with z on the + side of F exists; similarly for /-(<). 




Figure 3: Contour F = F^ U Tc U F^ with the orientation chosen. 

Also, we denote by Q-i and f2i the bounded components of C \ F containing —1 and 1, respectively, and 
by floo the unbounded component of C \ F (Fig. 
In what follows we write 

R{z) = {{z - c+){z - C-)y/\ zec\rc, (2.4) 

which is defined and analytic in the cut plane C \ Tc, such that R{z) ^ z a,s z oo. 

We also need the critical orthogonal trajectories of the quadratic differential H2.3|l . These are defined by 



where 

F^- = |z€C-: Im^' J^di = o' 
where the integration is along a path from C_ to z in \ Tc, and 



e C+ : Im / -^dt 



where the integration is along a path from C+ to z in C+ \ Tc- 

The typical structure of the orthogonal trajectories F''- is shown in Fig. ^ Three orthogonal trajectories 
emanate from both (^^ and ending at 1, —1 and oo, respectively (see the dotted lines in Fig. ^J. We 
denote by 7^*", j^i, 7^ the arcs of F^ that connect C+ with the points 1, —1, and 00, respectively; this is 
also the part of F"*^ in the upper half plane. The corresponding arcs in the lower half plane are denoted by 
7j~, jZi and 7^, so that 77 is the mirror image of 7^ in the real axis, for se{l,— l,cx)}. 

2.2 Weak convergence of zeros 

Then we define the absolutely continuous (a priori, complex) measure /i on F by 

, , , A + B + 2R+(z) , , , 

dK^) = — ^T^dz, zeT, (2.5) 
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where i?+ denotes the boundary value of R on the +-side of F. (Only on Tc there is a difference between 
the + and — boundary values.) The line element dz is taken according to the orientation of T. 

Lemma 2.2 The measure H2.5|l is positive and 

^(rL) = i + A>o, M(rc) - -1- A-B > 0, /i(rfl) = 1 + B > 0. (2.6) 

In particular we have that ^ is a probability measure on T. 
The importance of fi is shown in the following result. 

Theorem 2.3 Let (an) and (/3„) be two sequences such that a„/n — > A and /3„/n B where A and B 
satisfy Suppose that 

lim [dist(a„,Z)]^/" lim [dist(/3„, Z)]^/" = lim [dist(a„ + /3„, Z)]^/" = 1. (2.7) 

n — *oo n — >oc n — >oo 

Then, as n oo, the zeros of the Jacohi polynomial p^"'^"^^ accumulate on T and ji is the weak* limit of 
the sequence of normalized zero counting measures. 

The conditions (|2.7|l imply that a„, /3„, and a„ + /3„ are not too close to the integers. That such a condition 
is necessary is easily seen from the case when these numbers are in fact integers (cf. H1.5(l ^. 
To describe the general case, we need the function 

'^(^) = — ^ — (2-^) 

which is a multi-valued function. However, its real part is well-defined, and we see from the definition l|2.2(l 
that Y — {^z : Re 0(2;) = 0}. For every r we introduce the level set 

= {z e C : Re0(z) = r}. (2.9) 
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We note that by the selection of the branch in (|2.4I) , Re > in the unbounded region fioo and Re < in 
the two bounded regions Q±i. For r > 0, we have that Tr is a simple closed contour in floo, while for r < 0, 
we have that Tr consists of two simple closed contours, one contained in Qi and the other in ri_i. We define 
for r < 0, 

We choose the positive (^counterclockwise) orientation on each of the closed contours. All these contours 
are trajectories of the quadratic differential ()2.3|l . See Fig.|Slfor the trajectories. 




Figure 5: Some trajectories of the quadratic differential H2.3|l . or equivalently, some level sets F^, for the 
values A = -0.7 and B = -0.8. 

Finally, we introduce three numbers r^, r^, and ra+p and we assume that 

(2.10) 
(2.11) 
(2.12) 

It is easily seen that these numbers are non- negative and that the case ra = r/^ — ra+p = corresponds to 
Theorem 12.31 It is also easily seen that at least two of the numbers Tq, rp and ra+p should be equal, and 
if the third one is different, it will be greater than the other two. So we distinguish four cases in the next 
theorem. 

Theorem 2.4 Let (an) and (/?„) be two sequences such that an/n — *■ A and /3„/n — > B where A and B 
satisfy (|l.l(l . Suppose that there exist three numbers ra, rp, and Va+f} such that the limits H2.10|l . (I2.11|l . and 
(I2.12|l exist. Then the following hold. 

(a) If r^ ^ rp = ra+p, then the zeros of Pn"'^'^"'^ accumulate on F as n oo, and fi is the weak* limit of 
the normalized zero counting measures. 



lim [dist(a„, Z)]^'" = 

n — ^oo 

lim [dist(/3„,Z)]^/" = 

n—*oo 

lim [dist (an H- /^n,^)]^^'' = 
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(Jo) If^a < r„+/3, then the zeros of Pi"''''^"-^ accumulate on Tc UT,. where r ^ [r )/2 > and 

A + B + 2R+{z) ^ 

5 Tdz, z e Tc U 

is the weak* limit of the normalized zero counting measures. 

(c) If ra — ra+f3 < rfj, then the zeros of pj°'"-^"-'> accumulate on Tn U Tr,-i where r = (r^ — f/3)/2 < 0, 
and 

A + B + 2 Rjz) ^ 

is the weak* limit of the normalized zero counting measures. 

(d) If rp — r^+p < fa, then the zeros of p^""^^"^ accumulate on Tl U F^^+i where r = [rfj — rQ,)/2 < 0, 
and 

A + B + 2 Rjz) ^ 
zm z'^ — 1 

is the weak* limit of the normalized zero counting measures. 

Of course the statement of Theorem 12 .31 is a special case of part (a) of Theorem l2.4l We choose to mention 
Theorem 12.31 separately, since it represents the generic case. The statements of Theorem 12.41 are also valid 
along subsequences of N, if we assume existence of the limits (|2.1()|l - ()2.12|l as n ^ oo for n in a subsequence 
A of N. 

To illustrate the different phenomena that can happen we show some figures (Fig. Eland Fig. UJ. 




Figure 6: Zeros of Pi^q for a = —70 +10 ^, /3 = —80 + 10 ^, together with the set F corresponding to 
A = -0.7, B = -0.8. 



Remark 2.5 A general approach to the limiting zero behavior of polynomials satisfying a non-hermitian 
orthogonality property has been established in the works of Stahl and Gonchar-Rakhmanov ^Hj . These 
authors describe the limit distribution in terms of the equilibrium measure in an external field on a contour 
satisfying a symmetry property in C. Our contour F possesses this property, but the theorems of |30| and 
|18j are not applicable: an essential assumption in these papers is the connectedness of the complement to 
the contour. Nevertheless, the measure [i from 12.5|l is the above mentioned equilibrium measure on F in a 
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Figure 7: Zeros of P^^f' for a = -70 + lO^^o^ /3 = -80 + lO^^o (left), and a = -70 + 10"^ + lO-^°, 
(3 = -80 - 10"^ (right), together with the set T corresponding to A^ -0.7, B = -0.8. 

certain external field. Also the contours Tr have the symmetry property and the measures given in parts 
(b)-(d) of Theorem 12 . 41 are the equilibrium measures in the external fields on these contours. So Theorem l2.4l 
shows that in a certain sense the results of Gonchar-Rakhmanov-Stahl are also valid for Jacobi polynomials 
with varying negative parameters. It seems likely that similar results hold in more general situations. 

2.3 Strong asymptotics away from (± 

The weak convergence results of Theorems 12.31 and 12.41 follow from the strong asymptotic results that we 
obtain for the Jacobi polynomials. We state the result here for the sequence p^^"^-^")^ ^gg p(An,Bn) 
denote the corresponding monic Jacobi polynomial. 

Note that F and F"*- divide the complex plane into six domains, which we number from left to right by 
I, II, III, IV, V, and VI, as shown in Fig. |H1 

To state the asymptotic results we need to be specific about the branches of the functions that are 
involved. We already defined <j) in (|2.8|l as a multi-valued function. Now we specify that 

= J^' ^dt, z e C \ (Fc U 7i+ U U 7i) (2.13) 

where integration from C_ to z is along a curve in C \ (F^ U 7j*^ U 7!]^ U 7^ ) . Note that this definition prevents 
the curve from going around the cut Fc and also from going around one of the poles ±1. 
Near infinity, (p behaves like 

0(^)-^-^f^log^ + c + oQ) (2.14) 

for some constant c. This constant c will also appear in the asymptotic formulas below. 

In our formulas we will also see fractional powers (z - l)-'4«/2 and {z + l)--^"/2^ We wiU choose these to 
be defined and analytic in C \ (7^ U 7+ ) and C \ {jti U 7^), respectively, and to be positive for real z > 1. 

Finally, we define 
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and 



1 I f^-c-V^' fz-c+^'^'' 



2i \\z — J \ z — ^_ 

which are defined and analytic in C \ Tc- The fourth- roots are chosen so that they approach 1 as z — s- cx). 
We call these functions iVn and N12 since they will appear later as the corresponding entries of a matrix N. 
Now we have all the ingredients to state our main theorem. 

Theorem 2.6 Let A and B satisfy Hl.l|l . Then the monic Jacobi polynomials p^'^"-^") fiave the following 
asymptotic behavior as n —> 00. 

(a) For z in domains I and II, 

(b) For z in domain III, 



^Bnm smjAnn) ^ 



smiiA + B)mr) ' \n J J (2.16) 



j^^,^ sm{Bnn) , ^^1 



sm((A + BjuTT) ' ' " 



(c) Jbr z in domain IV, 

p(An,Bn)^^-^ = e-""(z - l)-^"/2(z + l)--B»/2 



-Amri 



Sin((.l + B)n7r)" '^"^'^ V V"/; (2-17) 

sin((A + i.)n.)^""'^^^^-(^)(l + K^)! 



_^gBn7r» sin(A717r) 
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(d) For z in domains V and VI, 

These asymptotic formulas hold uniformly for z in the indicated domains away from the branch points, 
uniformly for A and B in compact subsets of the region —1 < A < Q, —1 < i? < 0, —2 < A + B < —1, and 
for values of n such that {A + B)n is not an integer. 

Remark 2.7 One can verify that the asymptotic formulas (|2.15|l - (|2.18(l agree on the boundaries of the 
respective domains. 

Remark 2.8 The fact that the formulas H2.15|l - H2.18|l hold uniformly for A and B in compact subsets of the 
region given by (|l.l|l implies that we can allow varying values of A and B. In particular, we can consider two 
sequences («„) and (/3„) such that the limits H1.2|l exist and satisfy Hl.l|l . We then have asymptotic formulas 
for the Jacobi polynomials p^""''^"^ as in H2.15|l - H2.18|l with A replaced by An = OLnjn and B replaced by 
iinjn. Then we also have to realize that 0, c, A'^n and A^i2 are going to be n-dependent. Indeed, these 
quantities are defined using A and B, which here we have to replace by A^ and _B„; we chose to state the 
theorem for a„ = An and /3„ = Bn for the sake of brevity of notation. 

Remark 2.9 The expressions between brackets in the right hand-sides of (|2.15(I - H2.18(I contain two terms 
that correspond to the Liouville-Green approximation of two linearly independent solutions of the differential 
equation satisfied by the corresponding Jacobi polynomials (cf. pTl Ch. VI] ) . In different regions of the plane 
and depending on the parameters, these two terms are of comparable sizes (and then zeros of the polynomials 
arise), or one of them is dominating the other. If we assume that An, Bn, and {A + B)n are not close to 
integers, the expressions sin(^n7r)/ sin((A + B)mr) and sin(i3n7r)/ sin((A + B)m:) have moderate sizes (not 
too small, not too big). In that case the dominant term is determined by Re0. For z G fioo, we have 
Re0(z) > 0, and then (|2.15() and H2.18|l both reduce to 

^g-nc(^_ -^)-A«/2(^_^^)-i3„/2g«0(.)^^^(^) ^ Q [j^^ (2.19) 

for z in domains I and VI. 

For z e Oi U we have Re0(z) < 0, so that e-"'^(^) dominates e"'^^^) for large n. Then (|TT3|) - (|TTH|) 
reduce to 

p{Au,Bn)^^^ = _ l)-^"/2(z + l)-^"/2e-"^(-)7Vi2(z) (^ + (2.20) 

for z in domains II and III (that is, for z G r2_i), and to 

p(A„,B„)(^) ^ e^""-"^ siii(^n7r) _ + l)-^"/2g-n0(.)^^2(^) (i^o(-\\ (2.21) 

sm{{A + B)mT) \ \n J J 

for z in domains IV and V (that is, for z £ J^i). 

We emphasize that H2.19|l . H2.2()|l . and H2.21(l only hold if An, Bn, and {A + B)n are not close to integers. 
In general one has to use the compound asymptotic formulas H2.15(l - (|2.18|l . 

Remark 2.10 If An is an integer, then (|2.17|l and H2.18|l reduce to (|2.19(l for z in domains IV, V, and VI. 
Then we see the multiple zero at 2: = 1, not only because of the factor (z — 1)^"^"/^, but also because 

(/>(2) = --log(z - 1) + 0(1) asz^l 
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so that 

e"'^(^) = (z-l)-'^"/2(l + 0(z-l)) asz^l. (2.22) 

So we have a zero at z = 1 of multipHcity —An, as it should be. 

Similar remarks apply if Bn is an integer. In that case we have a zero at z = — 1 of multiplicity ~Bn. 

Remark 2.11 If An is not an integer, then does not have a zero at z = 1. This is in agreement 

with formulas H2.17|l and (|2.18|) since the zero at z = 1 due to the factor (z— l)-^"/2 ig compensated exactly 
by the singularity in e^"'^(^) at z = 1, see H2.22|) . 



2.4 Strong asymptotics near 



The asymptotic formulas (|2.15() and (|2.18|l are not valid near the branch points C- and C+- Near those points, 
there is an asymptotic formula involving Airy functions. We need the following particular combination of 
Airy functions, depending on A, B, and n, 



A{s;A,B,n) 



sin(A7T,7r) 



■ A w Ai(ws) + e 

sm((A + B)m:) 

1 cos((A + B)n'K) - exp((i? - A)m:i) 

2i sin((A + B)mT) 



-AnTTi 



sui{Bmr) 
sm{{A + B)mi) 

Ai(s) + lBi(s), 



uj"^ Ai(u;^s) 



(2.23) 
(2.24) 



where uj — e^'^^l'^ and Ai and Bi are the usual Airy functions 1 . Note that A[s\ A^B ,n) is defined for 
combinations of A, i?, and n that are such that (A + B^n is not an integer. 

Theorem 2.12 Let A and B satisfy (|1.1|) . Then there is a d > such that for every z with |z — C_| < S, 
the monic Jacobi polynomials p^(^"^^") fiay^ f/jg following asymptotic behavior as n oo: 



^1/6 



-j-An/2 

z-C- 

1/6 



(^ + 1) 



/(z)^ ' A{n"^f{z)-A,B,n)[l + o[^^ 



(2.25) 



-1/4 



/(z) ^{n'/^f{z);A,B,n)[l + 0[- 



with 



-1 2/3 



(^) 



(2.26) 



where the 2/3r(i root chosen is real and positive on 7^. The O -terms in (|2.25|l hold uniformly for |z — C_| < 6 
and for A and B in compact subsets of the region —1 < A < 0, —1 < B < 0, —2 < A + B < —1, and for 
values of n such that {A + B)n is not an integer. 

There is a similar asymptotic formula for the behavior near 

Remark 2.13 From the uniform asymptotics in Theorem 12.121 it is possible to establish a more precise 



behavior of the zeros of Pn 



(An,Bn) 



close to the branch points C±. In fact, the zeros of the function A defined 



in (|2.23|l - (|2.24() model the behavior of these zeros of p^^"-'^"-\ Yot: instance, in the generic case (12.71) both 
terms in H2.23|l or (|2.24|l have approximately the same size, and A has its zeros aligned along three curves 
emanating from and forming the same angle. The situation is different in cases (b)-(d) of Theorem 12.41 
For instance, if Va+p > ra = Va, then the first term in (|2.24|) dominates the second term. But if 7^ 
then one of the two terms in H2.23|l dominates the other. In these cases the zeros of A behave like zeros of 
the dominating Airy function and are aligned along a single curve emanating from 0. 
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3 Proof of Lemmas [O and 

Proof of Lemma 12.11 The quadratic differential H2.3|l has a simple zero at (± and a double pole at ±1 
and at oo. This determines the local structure of the trajectories as follows, see also |2H1 Chapter 8] or 
[SH Chapter III], 

(1) Three trajectories emanate from C± at equal angles. These are the critical trajectories. 

(2) Near ±1 the trajectories are simple closed contours. Here we use the fact that 

i^-C.)i^-C-)- ,o(^] as.^±l, 



(z2-l)2 (zTl)2 ' " V^Tl, 

with c±i < 0. 

(3) The trajectories near oo are also simple closed contours. This follows from the fact that in the expansion 

(^-(-)(^~^^)^2^^0(X] as. 



(z2_l)2 z2 

we have Coo = — 1 < 0. 

In the lower half-plane C~ there is only the simple zero at C_. The other points are regular points. 
This means that the three critical trajectories that emanate from C_ extend to the boundary of C~, cf. |28l 
Lemma 8.4]. Because of (2) and (3) and the fact that trajectories do not intersect, the critical trajectories 
do not tend to infinity, or come to ±1. So each critical trajectories exits the lower half-plane in a point from 
M \ {—1, 1} and these points are mutually distinct, say ^c, and with £_l < £,c < £,R- Because of the 
symmetry with respect to the real axis, 

Jc- ^'-1 

Hence, the three trajectories extend into the upper half-plane as their mirror images in M, and so they 
continue to This proves the existence of three arcs F^, Tc, and Tj^ contained in T and connecting (^±, 
where G for s e {L,C,R}. 

Next, we note that Fl U Tc is a closed contour consisting of trajectories. It follows from [281 Lemma 
8.3] that it has to surround a pole. Similarly Fc U Tji has to surround a pole. This can only happen if 

< -1 < < 1 < ^R- 

To complete the proof of the lemma, we need to establish that F consists only of F^, Fc, and Tfi and 
nothing more. We use that the function 

h{z)=ReJ^^ ^^dt, zeC\(FcU{-l,l}), (3.1) 

is single- valued and harmonic in C \ (Fc U { — 1, 1}). The path of integration in H3.H) is in C \ Fc. It is easy 
to see that 

lim h{z) = +00. 

z — ^oo 

Since /i = on F^, U — dflco, it follows by the maximum principle for harmonic functions that h(z) > 
for z e floa- Similarly, since 

lim h{z) — —oo, 

z — >oo 

and /i = on F^UFc = d^-\, and on Ffl;UFc = Sfii, we have that h{z^ < for z S ri±. Since V — {h — 0}, 
we get that F consists exactly of F^, Fc, and F^. This completes the proof of Lemma [2. II □ 
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Proof of Lemma 12.21 Recall that R{z) := — C+)(z — C_) denotes the single- valued branch in C \ Tc 
such that R{z) ^ z as z —^ oo. With this convention and taking into account H2.1() it is straightforward to 
check that 

O D 9/1 

i?(-l) = __ < 0, i?(l) = — - > 0. (3.2) 

^ ' A + B + 2 ^ ' A + B + 2 ^ ' 

From the definition of F it follows that dfi{z) is real- valued on F and does not change sign on each component 

ofF\{C-,C+}- 

Using the residue theorem, we have that 

/i(Fc U F^) = / dfi{t) = iA + B + 2) res (4^) = 
JFcurn z-1 \z - ly 

where we have used (|3.2() . Analogously, 

M(Fi U Fc) = / dfiit) ^iA + B + 2) res ( 



B. 



Finally, 



Hence, 



and 



U Tr) = [ dii[t) = {A + B + 2) res 



R{z) 



= A + B + 2. 



^i{T) = i {n{Tc U Tr) + ^(Fi U Fc) + ^^{Tl U F^^)) = 1, 



/x(Fi) = l-Ai(FcUFfi) = \ + A, 
M(rij) = l-Ai(FLUFc) = 1 + 5, 
^l{Tc) = l-M(rLUFfl) = -1-A-B, 

which proves (|2.()|) . Since each part has positive total /i-mass and /i does not change sign on each of the 
parts, we find that is a positive measure. This completes the proof. □ 

4 A Riemann-Hilbert problem for Jacobi polynomials 

Consider a closed path F„ encircling the points +1 and —1 first in the positive direction and then in the 
negative direction, as shown in Fig.|^ The point ^ G (—1, 1) is the begin and endpoint of F„. 




Figure 9: The universal curve F, 
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For a, /3 G C, denote 

w{z; a, (3) := (1 - z)"(l + = exp[a log(l - z) + /31og(l + z)]. 

This is a multi-valued function with branch points at oo and ±1. However, if we start with a value of 
'w{z;a,P) at a particular point of Tu, and extend the definition of it;(z;a,/3) continuously along r„, then 
we obtain a single-valued function w{z; a, (3) on r„ if we view as a contour on the Riemann surface for 
the function w{z;a,P). For definiteness, we assume that the "starting point" is ^ G (—1,1), and that the 
branch of w is such that w(^;a,/3) > 0. We prefer to view Tu as a subset of the complex plane. Then r„ 
has points of self-intersection, as shown in Fig. |5| At points of self-intersection the value of w{z; a, (3) is not 
well-defined. 

In (20) it was shown that for k € {0, 1, . . . , n}, we have 

_ 2 2n+a+/3+3 7ri(Q+/3) 

1 (2n + a + p + 2)1 (— n — ajl (—n — p) 

This shows that the Jacobi polynomials P,!"'^'' are orthogonal on the universal curve r„. The right-hand 
side of H4.1|l vanishes for k ^ n \i and only if either — 2rt — a — /3 — 2, orn-|-aorn-|-/3isa non- negative 
integer. In some of these cases the zero comes from integrating a single- valued and analytic function along a 
curve in the region of analyticity; other values of a and (3 correspond to the special cases mentioned before 
when there is a zero at ±1. It is shown in [201 that the orthogonality conditions 1)4. 1|1 characterize the Jacobi 
polynomial Pn"''^^ provided the parameters satisfy 

-n - a - /3 ^ N, and n + a ^ N, and n + (3 (^N. (4.2) 

Then pi"'''' is of de gree exactly n, and we will denote by Pn"'^^ the corresponding monic Jacobi polynomial. 

Based on the orthogonality (|4.1|) a Riemann-Hilbert problem is constructed in [^, whose solution is 
given in terms of Pn"'^^ with parameters satisfying (|4.2|l . 

Let r„ be a curve in C as described above with three points of self-intersection as in Fig.|51 We let r° be 
the curve without the points of self-intersection. Recall that the orientation of r„ (see also Fig. induces 
a -|- and — side in a neighborhood of r„, where the -I- side is on the left while traversing r„ according to its 
orientation and the — side is on the right. Again, we say that a function 1" on C \ r„ has a boundary value 
Y+{t) for t e r° if the limit of Y{z) a,s z ^ t with z on the -I- side of r„ exists; similarly for Y^{t). 

The Riemann-Hilbert problem for Jacobi polynomials is then as follows. We look for a 2 x 2 matrix 
valued function Y = : C \ r„ ^ C^^^ such that the following four conditions are satisfied: 

(a) Y is analytic on C \ r„. 

(b) Y has continuous boundary values on r°, denoted by F+ and such that 



(c) As z — > oo, 

Y{z)=(l + o' 







— n 



(d) Y{z) remains bounded as z — > t G r„ \ F". 



This Riemann-Hilbert problem is similar to the Riemann-Hilbert problem for orthogonal polynomials 
due to Fokas, Its, and Kitaev see also . Also the solution is similar. 
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Proposition 4.1 ([20 ) Assume that the parameters a, (3 satisfy H4.2|) . Then the above Riemann-Hilbert 
problem for Y has a unique solution, which is given by 



Y{z) 



Pl°'-f'>(t)w(t;a,p) 



dt\ 



dt 



ze c\r„, 



(4.3) 



r„ 



for some non-zero constant c„ 



The Riemann-Hilbert problem holds for any combination of parameters a and /? such that (|4.2|l is 
satisfied. Also the contour is rather arbitrary. It could be modified to any curve that is homotopic to it 
inC\{-l,l}. 



5 Transformations of the Riemann-Hilbert problem 

In this section we consider parameters A and B satisfying the inequalities p.l|l . We also assume that n e N 
is such that An, Bn and (A + B)n are non-integers. Throughout this section A, B, and n remain fixed. 

From Proposition l4. II we know that the Jacobi polynomial p^^"'-^") jg characterized as the (1, 1) entry of 
the solution of the Riemann-Hilbert problem for Y given in the previous section with a — An and f3 — Bn. 
In this section we apply the steepest descent method of Deift and Zhou to this Riemann-Hilbert problem in 
order to reduce it to a Riemann-Hilbert problem that is normalized at infinity and whose jump matrices are 
close to the identity. In the next section we derive the asymptotic results from this analysis. The Dcift/Zhou 
steepest descent method proceeds through a number of transformations of the original Riemann-Hilbert 
problem. 

5.1 Choice of contour 

In the first step of the analysis we have to pick the right contour. For A and B satisfying (|1.1() we have 
the contour F = F'^'^'^^ defined in H2.2|l . which according to Lemma l2.ll consists of three analytic arcs 
F = Fi U F(7 U Ffl. We modify F„ to a contour that follows F in such a way that every part of F is covered 
twice as shown in Fig. HJH 




Figure 10: Tautening F„ on the set F. 
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Passing from the Ricmann-Hilbcrt problem on r^, to the Riemann-Hilbert problem on F, we have that 
on each part of F two of the jumps are combined. The new jump matrices take the form 

1 w{ti;An,Bn)\ I'l -w{t2; An, Bn)\ _ f'l w{ti; An, Bn) - w{t2; An, Bn) 
1 j 1^0 1 ) ^ \S> 1 

where ti and <2 are points on the Riemann surface, both lying above t. The values of 'w{tj; An, Bn), j — 1, 2, 
differ from each other by a phase factor. To make this precise we specify a single- valued branch for the 
weight 

w{z; An, Bn) = (1 - z)^"(l + z)^" 

on F. Since F \ {C+} is simply connected, we can define a single valued branch on F \ {C+}, and we will do it 
in such a way that An, Bn) > 0, where ^ = f c is the intersection point of Fc with the interval (—1, 1). 

Then the jump on each of the contours F^, Fc, and can be calculated. The result is the following 
Riemann-Hilbert problem on F for a matrix valued function which we continue to call Y. The contour F has 
the orientation shown in Fig. |31 

(a) y is analytic on C \ F. 

(b) Y has continuous boundary values on F \ C+}i denoted by 1+ and Y-, such that 

Y+it)=Y^it)(^^ dswit;An,Bn)\^ for i e F, \ {C±}, se{L,C,R}, (5.1) 
with constants 

dL = e^'^-^™ ^^2^Am _ ^ dc^l- e2'^(^+^)", (Ir = 1 - e^'^'^™, (5.2) 
and we follow the convention about the branch of w{t; An, Bn) on F mentioned above. 

(c) As 2; — > oo, 



(d) Y(z) remains bounded as z (^±. 

Of course the solution to the above Riemann-Hilbert problem is similar to the solution H4.3|l to the earlier 
Riemann-Hilbert problem. In particular we still have 

rn(z) = Pr"^^")(z) (5.3) 

The constants d^, dc and d^ from (|5.2() will play an important role in what follows. These numbers are 
non-zero, exactly because of our assumption that An, Bn, and {A + B)n are non-integers. Observe that 

dL+dc = dn, (5.4) 

which is a relation that will be used a number of times. 

5.2 Auxiliary functions 

In order to make the first transformation of the Riemann-Hilbert problem we need some auxiliary functions. 
We already know from Lemma |2 . 21 that ^ defined in (|2.5|l is a positive measure on F such that 

^(Fl) = 1 + .4 > 0, ^l{Tc) = -l-A-B>{), ^l{TR)^\ + B>Q. (5.5) 
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Let g be the complex logarithmic potential of the measure ^, 



log(z — t)dfi{t). 



This is a multivalued function; however its derivative is single valued: 

A + B + 2 R{z) A/2 B/2 



9'{z) = 



dii{t) 



A + B + 2 R{z) 



z-1 z+1 
A/2 B/2 



z^ — 1 z — 1 z + 1 



for z € f^oo, 

for z e rj_i U fli. 



We define 



G{z) = exp g' 



it)dt], z e c \ r, 



where the path of integration lies entirely in C \ F except for the initial point C_ . From the fact that /i is a 
positive unit measure on F it follows that G is single- valued in each component of C \ F. Furthermore, G is 
analytic, G'(C-) = 1, and the following limit exists 



hm^.C-exp 



2— »00 Z 

We calculate the jumps of G. We have 



G+{z)G-{z) 



w{C-;A,B) 
w{z;A,B) 



and 



G+{z) 
G-{z) 



exp(-20+(z)), 



{g'{t)-l/t)dt 



for z e F, 



for z e F, 



where <j) is defined by H2.13(l . It will be useful to introduce also the related function 



A + B + 2 



m_ 

t2 



dt — 4'{z), for z G 

A+B+2 rC+ Bit) 



\(FcU7ZiU7rU7^). 



(5.6) 

(5.7) 
(5.8) 

(5.9) 



1 

To relate (/> with (p it is necessary to compute ^+^+^ ^iiijf^. This integral depends on the path from 



C- to C+- We can distinguish four paths, namely F/j, —Tc,+, —Tc,- and — F^. (Recall that Tc and F^ are 
oriented from ^+ to <^_. So we put a minus sign to indicate that the path is from to ('+.) We obtain 

ni^{Tfi) = ni{l + B) integral over F/j 

— 7ri^(Fc) = ni{l + A + B) integral over —Fc,+ 

Trifi{Tc) = — Trill + A + B) integral over —Fc,- 

— 7ri^(Fi) = — 7ri(l + A) integral over — F^ 



A + B + 2 



m 

t2- 1 



dt 



It follows that 










0+(z) = 4, 


z) 


+ 7ri(l - 








z) 


— 7ri(l - 


M) 




<^+(z) = 


'z) 


— 7ri(l - 






</._(z) = 


'z) 


— 7ri(l - 






0+(z) = 


'z) 


+ TTi{l - 






(/._(z) = 


'z) 


+ TTi{l - 







for z on 7+ , 
for z on 7^]^, 
for z on 7^. 



(5.10) 
(5.11) 
(5.12) 



Observe also that by construction both (p and 4> have negative real parts in the bounded components 17 _i 
and r^i of C \ F (where defined) and positive real part in (with the appropriate cuts). 
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5.3 First transformations Y U 

Now we introduce a new matrix valued function U by 

U{z) = k""^ w{C-;An, Bny^/^ Y{z) G(z)""'"3 w(C- ; An, Bny^"^ 



where '''3 = ( g Pauli matrix, and for any non-zero x, x"'^ = 



(5.13) 

Here G is the function 



introduced in (|5.2|l . and k is the limit defined in (|5.6|l . Then U satisfies the Riemann-Hilbert problem 

(a) U is analytic on C \ F. 

(b) U has continuous boundary values on F \ {C±} such that 



C/+(z) = C/_(z) 



exp(2n(/)+(z)) 




exp(— 2n(/)-|_(;z)) 



forzeF°, se{L,C,R}. 



(5.14) 



(c) U{z) = / + O (7) as z -> 00. 

(d) U{z) remains bounded as z — > C±. 

To obtain the jumps in H5.14|l we used the relations 1)5. 7|l and 1)5.8(1 . For the asymptotic behavior in (c) 
we used the limit H5.6|l . 

We use the following factorizations of the jump matrices in H5.14|l 



and 







-2nc 



dc 



1 







-e-2«^ 1 - 







in order to define the next transformation. 



dc 



for s ~ L, R, 



(5.15) 
(5.16) 



5.4 Second transformation U ^ T 

The trajectories F and the orthogonal trajectories F^ divide the complex plane into six domains, which we 
number from left to right as domains I, II, III, IV, V and VI, see Fig. |H1 We define T in each of these six 
domains separately. We put 



T 
f 

f = U 
f = U 
f = U 
f = U 



1 
1 

1 g-2n0 I 



1 



1 g2«0 I 

da 



1 g2n0 1 



1 e2"^ 1 



in domain I, 



in domain VI, 



-di 

-dR 
^ 



, -dL 



= C/l _i 



-rffi 



1 „2ri0 



0\ / 
1 



0\ / 



-dL 








in domain III, 
in domain IV. 



in domain II, 
in domain V, 



(5.17) 
(5.18) 
(5.19) 
(5.20) 
(5.21) 
(5.22) 
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Since 0(0) behaves like ^+^+'^ logz as z ^ 00, we have |e~2"'^("=)| ~ \z\ (^+B+2)n^ ^j^^^^ 



^g-2n0(z) ij^I + 0{l/z) asz^oo. 

Thus f(z) = / + 0(1/ z) as z -> 00. 

By definition, T is analytic in C \ (F U T-'-). However, we have arranged our transformation in a way that 
the jumps on Tl and Tn disappear (due to the factorization (|5.16() and the definition of T) so T is analytic 

inc\(rcur^). 

We compute the jumps on Tc U T-^ with the convention that these curves are oriented as shown in Fig. 
IHl Straightforward computations then show that 



V dLdn ^ J 



r+ = r_ I j_^_2„0_^ _ j_g_2«0_ on7^U7^, 

J ^ ^ J on 7+1 U 7-1, 

r+ = T_ ( J da + '^R^ \ on 7+ U 7r- 

Since is analytic across the curves , the jumps on these curves simplify to (we also use 1)5.4(1 ) 

~ ( I 0\ 

r+=T_(^^ '^^^ j on7ZiU7r. 

If we now express the jumps on the contours jj' in terms of cf), see (|5.10|l - (|5.12f) . they look as those on the 
lower half plane, but with (jj replaced by (fi: 

T+ = f^\ dn ^„~2n4> ? ) on 7+ 



dc ^-2n4> ij "n rooi 
dndi, / 

^ dhdR ^2na 



T+=T_(^ dc^'^ ) on7+iU7i+ 



Now with r such that 

2 dhdn 

we define T by 





(5.23) 



^-^0 rj^^O r-j- (^-24) 

The effect on the jump matrices is that the (1,2) entries are multiplied by and the (2,1) entries are 
multiplied by r^. So T satisfies the following Riemann-Hilbert problem: 

(a) T is analytic on C \ (Fc U T^). 
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(b) T has continuous boundary values on (Fc U F^) \ {C±} such that 



r+=r_(^_°^ Jj onFc, (5.25) 

T+ = r_ (^_J^2n^ on 7^, (5.26) 

T+=T-i^Q \ j on7ZiU7r, (5-27) 
1 

r+=r_fl -f"^) on7l,U7i+. (5.29) 



T^+=r-( J on7i, (5.28) 



1 



(c) T{z) = I + 0{l) asz^oo. 

(d) T(z) remains bounded as 2 ^ C±- 



The problem for T is by now relatively standard. However, compared with earlier works, the triangularity 

i 



of the jump matrices on the curves 7^ is reversed. The inverse transposed matrix T * satisfies the jumps 



T^'=TZ'{\ J) onFc, 



rj-l — t rj-i — t 



1 



on 7^, 



^ j on7i, 
T+'=TZ'[^^l~^ J) on7liU7i+, 
which are exactly of the form considered for example in [3 1221 • 
5.5 Outside parametrix 

The jump matrices in H5.26(l - (|5.29|l are close to the identity matrix if n is large. Therefore we expect that the 
main term in the asymptotic behavior of T is given by the solution N to the following model Riemann-Hilbert 
problem: 

(a) N is analytic in C \ Fc, 

(b) N+ - iV_ on Fc \ {C±}, 

(c) N{z) = 1 + 0{l/z) as z -> 00. 

In analogy with the condition (d) in the Riemann-Hilbert problem for T we would like to ask that N{z) 
remains bounded as z (l±. However, this would lead to a Riemann-Hilbert problem with no solution. 
Instead we allow for moderate singularities of at 
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(d) N{z)^0{\z~(:±\-^/^) asz-.C±. 
The solution to this problem is given by 



au)^a(.r^ au)^(.r^ (5-30) 

2i 2 / 

with 

being analytic and single- valued in C\Tc, such that a{z) ^ 1 as z oo, see [H1|3|S1E2]. In [52] also an 
alternative expression for TV has been established in terms of R{z) :— ^y (z — (^){z — (^): 

1/2 /i_p',.,\l/2N 



^W-l 1,1/2 .1.' /i/2 I- (5-31) 



I--R'(Z) \ ' l + fl/(2) 



5.6 Local parametrices 

Near the branch points (± we construct local parametrices in the same way as done by Deift et al |3 13 E] , 
see also 122) • In a neighborhood Us = {z C : |z — C_| < (5} of ^_ we construct a 2 x 2 matrix 
valued P that is analytic in Us \ (Tc U •yZi U 7f U 7^), satisfies the same jump conditions as T does on 
Us n (Tc U U 7f U 7oq) and that matches with N on the boundary Cs of Us up to order 1/n. 




Figure 11: Conformal mapping /. 



We need the function 



/(^) 



2/3 



(5.32) 



where the 2/3rd root is chosen which is real and positive on 7^. This is a conformal map from Us onto a 
neighborhood of provided 5 > is small enough. We note that 7^^ is mapped to the positive real axis, Tc 
to (a part of) the negative real axis. Recall that is real and negative on 7^" and 7Z1 and we see that 7]~ 
is mapped to argw — 2tt/3 and 7Z1 to argw — —2tt/3 (Fis. Ill|l . Then the Ricmann-Hilbert problem for P 
is solved by (cf. |12]) 

P{z) = [i;(z)*(n2/V(2))e"'^*'^^'''l ' , (5.33) 
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where 



E{z) 



1 -1 

— i —i 



a{z) 



(5.34) 



and is built out of the Airy function Ai and its derivative Ai' as follows 

for < args < 27r/3, 



Ai(s) Aiiuj^s) 

Ai'(s) uj^ Ai' (uj^s) 

Ai(s) Ai(w2s) 

Ai'(s) uj^Ai'iuj^s) 

Ai(s) -Lu'^ Ai{ujs) 
Ai'(s) - Ai' {lus) 

Ai(s) -uj'^ Ai{ujs) 
Ai'(s) -Ai'(ws) 



e 6 



e 6 



e 6 



1 

-1 1 

1 
1 1 



for 27r/3 < args < tt, 
for — TT < args < — 27r/3, 
for — 27r/3 < args < 0, 



(5.35) 



with uj = e2^^/3. 

Note that we take the inverse transpose in H5.33|l . which is absent in the construction in This is 

of course due to the fact that the Riemann-Hilbert problem for T~* is comparable to the Riemann-Hilbert 
problem found in [221, see the remark at the end of subsection 5.4. 

A similar construction yields a parametrix P in a neighborhood Us — {z : |z — C+| < S}. 
5.7 Final transformation T \—>- S 




Figure 12: Contour Ts for the Riemann-Hilbert problem of S. 

The final transformation T t-^ S is 

S ^ TN^^ outside the disks Us and Us, 
S = TP-^ inside the disk Us, 
S = TP-I inside the disk tjs- 



(5.36) 
(5.37) 
(5.38) 
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Then by construction, S has jumps on the circles Cs — dUs and Cs = dUs as well as on Tc U F^. Since 
the jumps of T and TV on Tq agree, we have that S is analytic across the part of Tc outside the disks Us 
and Us- Similarly, the jumps of T and P agree inside the disk Us, and the jumps of T and P agree inside 
the disk Us, so that S is analytic in Us and Us with possible isolated singularities at C±. However it follows 
from the behavior of T and N near C± that the singularities are removable. Thus S solves the following 
Riemann-Hilbert problem. 

(a) S is analytic on C \ Ts, where Ts consists of the circles Cs and Cs, and of the parts of 7_i, 71 and 700 
outside the disks, see Fig. 1121 

(b) S has continuous boundary values on F5 such that 

S+^S 
S+^S 

s+^s 

s+ = s 

s+ = s 
s+ = s 

(c) S{z) = / + (i) as z ^ 00. 



-^(^_g-2„0 ^jN-' on J- \Us, (5.39) 

-N L \ j on (7Z1 U 7r) \ Us, (5.40) 

-^(_gi„0 ^^N-' oni+\Us, (5.41) 

-N (^J \ j iV-i on (7^l U 7i+) \ Us, (5.42) 

_PA^^i on Cs, (5.43) 

_PA^^i on Ca. (5.44) 



6 Asymptotics: Proofs of the theorems 
6.1 Asymptotics of 5" 

The analysis in the last section is done for fixed values of A, B, and n. All the transformations are exact for 
finite n. It is now our aim to let n — s- 00 and control the jump matrices in the Riemann-Hilbert problem for 
S. We want to do it in a way which is valid locally uniformly for parameters A and B satisfying (|l.ll) . 

Then first of all we should study the dependence of the contour Fg on the parameters A and B. Note 
that Ts does not depend on n, but it does depend on A and B. In fact, we have that F^ is completely 
determined by A and B, while the radius S of the circles around (± is only restricted by the requirement 
that the mapping / from (|5.32|l is a conformal mapping on Us- With that in mind, it is clear that we may 
assume that the curve Ts depends on A and i? in a continuous way. 

Now we can see what happens with the jump matrices in (|5.39|l - (|5.44|l as n ^ 00. On 7^ \ Us we have 
that Tlecf) is strictly positive. Hence the jump matrix in (|5.39() is / + 0(e~™) as n — > 00, uniformly on 
7oo \ Us- This estimate is also valid uniformly for A and B in compact subsets of the set 

{{A, B)\-l<A<0, -KB <0, -2<A + B< -1}. (6.1) 

Similarly, the jump matrices in H5.40|l - H5.42() are /-l-0(e^'^") as n — s- 00, uniformly on the respective contours 
and uniformly for A and B in compact subsets of (j6.1(l . 

For H5.43|l and H5.44|l we make use of the matching conditions 

P{z) = \I + o(-]] N{z) uniformly for z G Cs. (6.2) 
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and 

P{z) = (l + o(-] ] N{z) uniformly for z e Cs- (6.3) 



So that the jump matrices in (|5.43l) and (|5.44() are / + 0{l/n) as n — > oo, uniformly on the two circles. A 
closer analysis also reveals that the O-terms in (|6.2ll and H6.3|l are valid uniformly for A and B in compact 
subsets of (|6.1|) . 

So the conclusion is that all jumps in (|5.39(l - (|5.44|l arc / + 0{l/n) uniformly for z on Ts, and uniformly 
for A and B in compact subsets of (|6.1|l . Then arguments such as in |S1[71|S] lead to the following conclusion. 



Proposition 6.1 We have that 

5(z) = / + 0(^) (6.4) 

uniformly for z (z C\Ts and uniformly for A and B in compact subsets of the set H6.1|l . 

The estimate (|6.4|l is the basic asymptotic result. Unravelling the sequence of transformations Y U t-^ 
T 1-^ r I— > S", we obtain asymptotic formulas for Y in any region of the complex plane. In this way we obtain 
the asymptotic formulas for P„ = Yu. 

6.2 Proof of Theorem ITf>l 

Proof of Theorem 12.61 Suppose A and B satisfy Hl.l|l and let n G N such that An, Bn, and {A + B)n 
are non- integers. Then we have Pi^"'^")(z) = Yii(z) by (Oil . 

For z in domain I away from the branch points, we get by using (|5.13|l . H5.17|l . I|5.24|l and H5.36|l . 

Yr,{z) = (^)"c/n(z) 

fG{z)\ frf, 1 „-2n<i(2 



. , Tn(z)-— e-^"*^^JTi2(z) 

\ K J \ dL 

- (^^^)"(e"^^^^Tn(.)-|e-"^(^)T.(.)^ 

= (^e"'^W(57V)n(z)-^e-"'^("H^^)i2(^)) • (6.5) 

Since S = I + 0{-^) and since the entries of N are bounded and bounded away from zero away from the 
branch points, we have that 

(5iV)ii = A^ii(/ + 0(l/n)) and (S'A^)i2 = Afi2(/ + 0(l/n)). (6.6) 

Next we recall that for z in domain I, 

G'{z) _ A + B + 2 R{z) A/2 B/2 

'g{z) -^^""'^ 2 Y - ^ 

so that 

A B 

log G{z) ~ (j){z) — — \og{z ^ 1) ^ log(z + 1) + const. 



Thus there is a constant c such that 

G(z)e-"^W 



e^^iz - \Y^I'^{z + 1)^-^/2 for z in domain I. (6.7) 
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Since 1^11(2:) is a monic polynomial of degree n, we find by letting 2 ^ cx) in H6.5|l and using (|2.14|) and (|6.7|l . 
that c should be as defined in H2.14() . 

Plugging H6.6() . H6.7() and the formulas (|5.2I) for d^j and into formula H6.5() we obtain H2.15|l for z in 
domain I. 

For z in domain II away from the branch points, we find in the same way 



) = (^^^^) (e-'^^'HSN)n{z)-^e--^(^\SN)M^ (6.8) 



Since G-\- — G—C ^"^ on Tl, see (|5.8|l . we have that Ge'^ is the analytic continuation of Ge into domain 
II. So we have by H6.7(l 

^^^^^^ =e-^(z-l)-^/2(^^^)B/2 for z in domain II. (6.9) 

K 

Then using (|5.2(l . H6.6(l . and H6.9|l in H6.8(l . we obtain H2.15|l for z in domain II. 
For z in domain III away from the branch points, we obtain 

^G(z)e^W V f_dL n^^,^, , , . _ dn n^i. 



yii(2) = [^^^ j (^-^e"^^-^(57V)n(z)-ge-"^^^^(5iV)i2(z)j . (6.10) 

Again we use (15. 2|) . I|6.6|) . and (|6.9|) to obtain (|2.16l) from (|6.10|) for z in domain III. 

The proofs of the formulas (|2.17(l and H2.18|l for z in the other domains IV, V, and VI are the same. 

We have derived the formulas H2.15|l - H2.18|) under the assumption that n is such that An, Bn, and 
{A + B)n are non-integers. Since the formulas hold uniformly in A and B in compact subsets of (|6.1|) and 
p{An,Bn) ^jgpgj^^g coutinuously on A and i?, they continue to hold if An or Bn is an integer. However, we 
cannot allow [A + B)n to be an integer, since then there is a reduction in the degree of P^^"'^") and we 
cannot normalize the Jacobi polynomial to be monic. 

This completes the proof of Theorcm l2.6l □ 

6.3 Proof of Theorem VTV^ 

Proof of Theorem I2.12L Suppose A and B satisfy (11.1(1 and let n e N such that An, Bn, and {A + B)n 
are non- integers. Then we have Pi^"'^"^(2) = Yii{z) by (lOll . 

Let z £ C/5 be in domain VI. Following the transformations ((5.13|l . ((5.2211 . ((5.24(1 . we see that 



Y,,{z)\ _ (G{z)e 



4>{z) \ " / pn<Piz) 



r(z) 



\ « / \ dc'^ 

where * denotes an unimportant unspecified entry. For z in domain VI, we have ((6.7(1 so that 



l^^iij^)^ = e-"^(z - l)-^"/2(z + 1)-^"/2t(z) 



d-L „-n(t>(z) 
dc 



Since z belongs to Us, we have T{z) = S{z)P{z) by (j^^Tzjl . By fHl?^ we have P(z) = i;-*(z)*^*(s)e-"'^(^)'"3 
where s = n^^^ f{z). Thus 

(^^"J^^^ = e-"=(z - ly^'^'^z + l)-^"/25(z)£;-*(z)*^*(s) ^ . (6.11) 
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From H5.34|l we see that 



Furthermore, we have < args < n/3 for s — n^^^f(z) since z is in domain VI, so that we use the formula 
(|5.35|) to evaluate vE'~*(s). Taking into account that 

Ai(s) Ai{uj^s) \ ^ „/6 



^""^ ' Ai'(s) a.2Ai'(L^2^) j - 2^^ 



we have from (|6.12|) and H5.35|l 

"^""^'^ 1 lJ\y^J \-Ai{uj^s) Ai(s) 



Plugging this into l|fi.ll|l we get 



M'A=e--%z-l)-^-/'{z+l)-''-/'V^S{z)(~: l)(^y^ /-(-^Ai(.)+^2Ai(c.^.))' 
V * / ' V V^i/V V -^Ai(s)+c^2Ai(^2^) 

where the prime denotes the derivative with respect to s. 
Comparing with H5.2|l and 1)2.24(1 we see that 

A(s) = A(s; A, B, n) = -^AUs) + uj^ AUcu'^s). 

dc 

Thus 



(6.13) 



We derived the formula H6.13|l for z Us in the domain VI. Similar calculations for z G Us in the other 
domains give the same result, so ((6.1311 is valid in the full neighborhood Us of Now it remains to recall 
that Yii = p^^">^") and that S = I + 0{l/n) as n ^ CX3 in order to obtain 1(2.25(1 . 

We have derived 1(2.25(1 under the assumption that An and Bn are non- integers. Since the formula holds 
uniformly for A and B in compact subsets of ((6.1(1 and p^^^"-^") depends continuously on A and B, they 
continue to hold if An or Bn is an integer. 

This completes the proof of Theorem 12. 121 □ 

6.4 Proof of Theorem HH 

Proof of Theorem 12. 4L Conclusions of Theorem 12.41 are based upon the strong asymptotics obtained in 
Theorem El 

We let (a„) and (/?„) be two sequences such that a/n — > A and /3„/n _B where A and B satisfy the 
inequalities ((1.1(1 and we assume that the limits 1(2.10(1 - 1(2.1211 exist. 

Taking into account that formula 1(2.15(1 is uniform in A, B, for z in domains I and II we have 



sin(M .-n,i.)^^^^^J,^o(^' 



sin((Q!„ + /3„)7r) 
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Since the first factors in the right hand side have no zeros in domains I-II, z is a zero 

sin(/3„7r) Ni2{z) 



sin((a„ + /3„)7r) A''ii(z) 



1 + 



(6.14) 



But N12/N11 is uniformly bounded and uniformly bounded away from zero, if we stay away from the branch 
points C±. Thus, taking the absolute values in (|6.14|) . we see that the zeros in domains 1-11 away from the 
branch points must satisfy 



2Re (h(z) 



■log 



sin(/3„7r) 



sin((a„ + /3«)7r) 



O 



Analogously, the zeros of p^""'''") in the other domains III-VI away from the branch points satisfy 



2 Re = - log 

n 



sin(^„7r) 



2Rc <j>{z) = - log 



sin(a„7r) 
sin(Q;„7r) 



O 



2 Re (/>(z) = - log 



sin(/3„7r) 
sin(a„7r) 



oil 



sin((a„ + /3„)7r) 



oil 



for z in domain III, 



for z in domain IV, 



for z in domains V-Vl. 



Furthermore, for any sequence of real numbers (Kn), 



lim— log |sin(7rK„)| = lim — log |dist(K„, Z) 

n n n 72 



whenever either one of these limits exists. Thus, under the assumptions of the theorem, the zeros of Pn 
away from the branch points must satisfy 



(a„,/3„) 



2Re0(z) =r + 



fa+[3 — f[3, for z in domains I and 11, 

fn — ff3, for z in domain III, 

ff3 — fa, for z in domain IV, 

Sa+I3 — fa, for z in domains V and VI. 



(6.15) 



From the definition of the constants H2.1()|l - I|2.12|l it follows that the "generic" case is 

fa ^f/3 ^ fa+P ■ 



(6.16) 



Recall that Re 0(z) > in domains I and VI, and Re (f){z) < in domains 11-V. Hence by H6.15|l . if z e C\ F, 
then p/i""'^"-'(z) 7^ for sufficiently large n, which proves that the zeros can accumulate only on F. 
Next assume we are in case (b) of the Theorem 12 .41 that is, 

ra+f3 > ra=rp ■ 

By (|6.15() . the zeros cannot accumulate in domains II, 111, IV and V, nor on F^ U F/f. In domains 1 and VI 
they must satisfy 



Re0(z) =r + I - ) , r = 



r^+p - f/3 



fa+f3 - fa 



showing that they must accumulate on the curve F^ defined in H2.9|l . Hence, in this case the accumulation 
set belongs to Fc U F,.. 
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The rest of the cases is analyzed in the same fashion. 

Once we have estabhshcd that the zeros accumulate along curves in the complex plane, it remains to find 
the asymptotic zero distribution. We can use the differential equation (see e.g. [321 §4.22]) 

(1 - z^) y';^{z) + [/3„ - a„ - (a„ + /3„ + 2)z] ^^(z) + n{n + a„ + /?„ + 1) y„(z) = 

satisfied by j/„ — p^"-'P"\ Rewriting this equation in terms of /i„ = y'nl{''^yri) reduce it to the Riccati 
form 

(1 - z ) - /i„(z) + /i„(z) + h^{z) + = . (6.17) 



Let Vn denote the normalized zero counting measures of pj°^"-'l^"\ Using the week compactness of the 
sequence (I'n) we may take a subsequence A C N such that converge along A to a certain unit measure v 
in the weak*-topology. By the discussion above, u is supported on a finite union of analytic arcs or curves 
(level sets F^), and every compact subset of C\supp(;/) contains no zeros of pj^"'^"^ for n sufficiently large. 
Hence, 

Jy z — t z — t 

locally uniformly in C \ supp(j^). Taking limits in 16.17|l we obtain that h satisfies the quadratic equation 

(1 - z^) h^{z) + [B- A-{A + B)z] h{z) + A + B + 1^0, 

so that 

r diy{t) A + B + 2 Riz) I f A B \ 

- ' ' 'zeC \ supp(j/) , 



:-t 2 z2-l 2Vz-l z+1 

with R defined in l|2.4|l and C,± given in l|2.1|l . By Sokhotsky-Plemelj 's formulas, on every arc of supp(i^), 

d.iz) = ^±^^dz. (6.18) 

2711 z^ — 1 

Assume that (16.16(1 holds, so that supp(j^) C F. By 1)6. 18() . v' — /i' almost everywhere on supp(;/). Taking 
into account Lemma 12.21 and that v \s a, unit measure it follows that v — fj.. 

If Ta+p > fa = then supp(^) C Fc U F^, r = {ra+p — fa)/"^ > 0. Again taking into account Lemma 
E21and the normalization of v it follows that supp(i/) = Fc U F^. 

The remaining cases are analyzed analogously. This completes the proof of Theorem l2.4l □ 
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